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Abstract 

We consider the light-cone (LC) gauge and LC quantization of the dimensional 
reduction of super Yang Mills theory from four to two dimensions. After integrat- 
ing out all unphysical degrees of freedom, the non-local LC Hamiltonian exhibits 
an explicit N = (2, 2) supersymmetry. A further SUSY-preserving compactification 
of LC-space on a torus of radius R, allows for a large-A^ numerical study where 
the smooth large-i? limit of physical quantities can be checked. As a first step, 
we consider a simple, yet quite rich, "Coulomb approximation" that maintains an 
N = (1, 1) subgroup of the original supersymmetry and leads to a non-trivial gener- 
alization of 't Hooft's model with an arbitrary -but conserved- number of partons. 
We compute numerically the eigenvalues and eigenvectors both in momentum and 
in position space. Our results, so far limited to the sectors with 2, 3 and 4 partons, 
directly and quantitatively confirm a simple physical picture in terms of a string-like 
interaction with the expected tension among pairs of nearest-neighbours along the 
single-trace characterizing the large- A'^ limit. Although broken by our approxima- 
tion, traces of the full N = (2, 2) supersymmetry are still visible in the low-lying 
spectrum. 
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1 Introduction 



Supersymmetric (SUSY) gauge theories give us a useful playground where we can test and 
try to understand non-perturbative features of gauge-theories. On the other hand large- iV 
expansions of various kinds [H E] provide further simplifications in the non-perturbative 
dynamics to the extent that suitable combinations of SUSY and large- often result in 
models that can be (almost) fully understood analytically and/or numerically. A well 
known example of such a powerful mix is, of course, the AdS / CFT correspondence [3l H] . 

Another application of the same set of ideas is the so-called planar equivalence O |6] 
between gauge theories with Dirac fermions in the antisymmetric (or symmetric) 2-index 
representation of SU{N) and theories with Majorana fermions in the adjoint. Such a 
correspondence (which holds at large volume for a "common sector" in the two theories) 
could be physically relevant if iV = 3 is "large-enough" since the former theories become, 
in that case, QCD-like. On the other hand, the adjoint-fermion theories become, in 
some cases, supersymmetric hence allowing to make predictions such as that of the quark 
condensate in QCD from the known value of the gluino condensate in A/" = 1 super- Yang- 
Mills theory [71 E]. 

Furthermore, SU{N) gauge theories with Nj > 1 adjoint fermions appear to present 
another peculiar large- A^ virtue, volume independence O |T0], according to which their 
properties at infinite volume are still captured if we consider them on a partially com- 
pactified space M.^ x with periodic boundary conditions both for the bosons and for 
the fermions. As the center symmetry appears to be unbroken in this case even at small 
volume pjj, we can send the volume of the to zero and yet recover, at large- iV, the 
large-volume physics (unlike the case of the well known Eguchi-Kawai reduction [12] for 
fundamental fermions). Combining planar equivalence and volume independence we may 
dream of calculating part of the spectrum of QCD (modulo 1/A^-corrections) by a small- 
volume calculation in a gauge theory with Majorana adjoint fermions. 

Having in mind these motivations, we combine in this work large- reduction and 
SUSY ideas by considering the well-known Af = 1 supersymmetric SU{N) Yang-Mills 
(SYM) theory compactified to two dimensions. However, since the precise connection of 
this model with four-dimensional SYM theory is not yet clear [9], we should consider for 
the time being this model for its own sake as a two-dimensional (and, to our knowledge, 
so far unsolved) model. 

We will be using a LC gauge and LC quantization [13], which has the advantage that 
all the unphysical degrees of freedom can be integrated out at the price of introducing 
Dirac brackets, a common procedure for constrained systems. The final outcome is a 
2-dimensional = (2, 2) supersymmetric gauge theory on the light-cone. It appears to 
have an exact SUSY vacuum and a well-defined, normal-ordered, positive-semidefinite 
LC Hamiltonian that can be written as an anticommutator of supercharges. At large 

this Hamiltonian acts on color-singlet physical states that can be written in terms of 
single-trace operators: 

Tr(oi(A;i)t...o„(A;„)t)|0) 
where the different o\ will represent either bosons or fermions with different light-cone 
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momenta ki > 0. 

The full Hamiltonian is still quite involved. In this paper, after having identified some 
"leading" terms which present potential Coulomb-like infrared singularities, we will define 
a natural Coulomb approximation to the full Hamiltonian. We show that, in this rather 
drastic approximation, the model becomes a partially (i.e. A/" = (1, 1)) supersymmetric 
generalization [T3] of 't Hooft's model in two dimensions [TS] with non-interacting sectors 
characterized by the number and species of partons they contain. 

Some analytic results will be presented but most of the calculations will be numerical. 
This is done by discretizing the LC momenta |16j| ki = ne, thereby reducing the problem 
to the diagonalization of a large M x M hermitean matrix. We then prove that our 
discretization indeed converges to a continuum- limit once we take the e — > {M oo) 
limit. 

The rest of the paper is organized as follows: in Section 2 we sketch the derivation of 
the LC two-dimensional theory starting from the SYM Lagrangian in four dimensions. We 
give the explicit expressions for the SUSY charges, light-cone momentum and Hamiltonian 
that satisfy the usual SUSY algebra. In Section 3 we introduce our Coulomb approxima- 
tion and show how its apparent infrared divergences neatly cancel for colour-singlet states 
leaving behind the effects of a confining linear potential. In Section 4 we introduce a LC 
compactification of the two-dimensional theory as a devise to work with discrete values of 
the light-cone momenta. In Section 5 we give our analytic and numerical results for the 
2, 3 and 4 partons sectors of the Coulomb Hamiltonian and discuss the left-over traces of 
the original supersymmetry. Finally, in Section 6, we give some conclusions and a short 
outlook. 



2 Light-cone SYM4 and its dimensional reduction 



Our starting point is 4-dimensional J\f = 1 super Yang- Mills (SYM) theory as defined by 
the Lagrangian [T7]: 



C = --F'^^F''^"' + iX'^a^'D^X' = Tr 



The gauge group is SU{N) and we have set to zero a possible 6'-term since it can be 
rotated away in the absence of a gluino mass (i.e. as long as we do not break SUSY). We 
shall use the following conventions: 

d^A^ - d^A^ + ig[A^, A^] , D^X = d^X + ig[A^, X] . 



A — Aarpa p 



pa pa 
/II/ fund 



■^S- and fix the LC gauge A_ =0. As a 



We now introduce LC coordinates 
result we can rewrite ([1]) as: 

£ = Tr [{d^A+Y + 2{d^A{d+A^ - DA+) + d^A^{d+A - DA+))+ 



(DA - DA^ - ig[A\ A]f + 2V2iX^ 



-D 



X 



(2) 
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where we have introduced A = {Ai + iA2)/V2, At = {Ai - iA2)/V^, I^i+i2 = V^D = 
V2{d + ig[A, x]) and its hermitean conjugate -Di_i2 = V^D = \/2{d + ig[A^ , x]). 

The fields A^ and Ai become non-dynamical (in the sense that their equations of 
motion do not involve derivatives with respect to the lightcone time x~^) and can be 
integrated out. Furthermore, we can consider the reduction of the theory to = 2 by 
discarding all dependencies on x"*" and hence by setting to zero d and d. As a result the 
Lagrangian under consideration becomes: 



CL =Tr 



J+ + 2{d^Ad+A^ + d+Ad^A^) - g^[A, A^]^+ 



+ 2v^z(A;a+A2 - g'[A, X;]^[A\ A^]) 



(3) 



where the subscript io means that we have integrated out the non-dynamical fields and 
the superscript r tells us that this is now a theory in 2-dimensions x^. The (reduced) 
current turns out to be: 

J+ = tg[A\ d.A] + ig[A, d_A^] + ^2 ^^{A;, As} . (4) 

To obtain the Hamiltonian from (|3]) we have simply to perform a Legendre transform: 



iZI=Tr 



- J+ + / [A, A^f + 2V2ig^[A, \l]^[A\ \, 



(5) 



We have to keep in mind, however, that the right commutation relations are those 
obtained from Dirac -rather than Poisson- brackets: only then the constraints implicitly 
used will be preserved by the LC-time evolution generated by Hl^. This leads to the 
following quantization of the dynamical fields A and A2 in terms of creation/annihilation 
operators with a Fock vacuum |O)0: 



A"(0,x~ 
A^(0,x- 



dk+ 



/2^V2AF 
dk+ 



a\k^ 



+ b^''{k+)e+"'^' 



r(A;^ 



+ g^\k 



10 21/4^27? 

[a%k+),a^\k'+)] =r,6t^] = {r,^} = {g^,g^'} = 5{k+ - k'+)5^' 



(6) 



and of course all other (anti-) commutators vanish. Note that in order to obtain precisely 
the Dirac commutation relations we have to take our LC momenta k'^ to be positive. 

The theory thus obtained can be seen as the M = (2, 2) two-dimensional supersym- 
metric theory that follows from the dimensional reduction of A/" = 1 super Yang-Mills on 



^We have assumed a vanishing VEV for the scalar field although, even in Z? = 2, we expect the 
classical moduli space to be preserved since tunneling should be suppressed in the large- iV limit. If so 
the spectrum will depend on the particular point chosen in moduli space. One of us (GV) is grateful to 
A. Armoni and A. Schwimmer for a discussion concerning this point. 
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X T^. Note, however, that due to the non-susy invariance of the gauge fixing the trans- 
formations generated by the supercharges through Dirac brackets will not be the usual 
one but will have to be supplemented by an additional gauge transformation needed to 
restore the gauge constraint y4_ = 0. This is already true before compactification (where 
one gets M = 1 SUSY in D = 4) and keeps working after compactifying on X T i.e. 
if we keep only the zero-modes (hence setting d/dx-^ = 0). A straightforward calculation 
leads to the following form for the supersymmetric charges, momentum and Hamiltonian 
operators: 



dx-2V2d-A^\l'' 
-l={g2, Q2] = J dx- (29_At"a_yl" + iy2A*/9_A^) 



'k2 I ' 



(7) 



dk k 



Qi = J dx-{-2V2gr^''A^''d-A^ + 2igr^'\l\l^)^ 



\ *c 
^2 



2V4 



dkidk2dk^ S{.)f 



abc 



ksVh 



a\ b 



k29k3 ^ki^k2fl 



ks 



k^.x/ko 



{ks + ^2) / at 6 c ra fO\ n 
-T ^, , [yki9k29ki JkiJk2Ji 



2k^k2 
1 



ks 



{k2 - ks) 
k3k2 



ks 



fkxfk29ka + 9k[fk29k: 



at M c 



2V2 



{Qi,Qi}= / dx-Hl 



(9) 
(10) 



where 6{.) is a shorthand notation for 6{±ki ± ...) where the sign plus is for creation 
operator while the minus sign otherwise. We did not report here the explicit momentum 
space version for P+ since, as the reader may guess from the rather involved expression for 
Qi, its expression is quite long and not particularly illuminating. A lengthy calculation 
shows that the above charges satisfy the J\f = (2, 2) SUSY algebra: 



{Qi,Qi} =2V2P+ = 2V2HI 
{Q2,Q2}=2V2P-, 



fill 



with all other anticommutators vanishing identically. 

As we will see in the following sections one can compactify also the x~ direction and 
replace all the integrals dk with sums over positive integers ^^1, the key point is 
precisely that all LC momenta have to be greater than zero, so we can actually rewrite 
Jq°° dk = dk 9{k). We note, incidentally that in [18] it was claimed that the discretized 
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version of the supercharges does not satisfy the susy algebra. We claim instead that 
everything works fine provided an appropriate care is used in performing the discretization 
and in defining the (anti) commutators; particular attention is needed when, by momentum 
conservation, an intermediate parton gets a vanishing momentum: by adopting a careful 
prescription for the ensuing ^(0) we can fully maintain supersymmetry. 

2.1 General properties of the reduced theory 

Let us discuss some exact features of the complete SUSY charges and Hamiltonian. The 
Hamiltonian conserves parton number parity ( — 1)^ (actually changes p by 0, ±2) and 
therefore splits into two sectors with even and odd number of partons. On the other hand 
the charges Q2, Q\ preserve p while Qi,Q\ change p by ±1. 

Besides P_ there is another conserved quantity, J^, which is what remains of the 
original helicity of the 4-D theory. The four partons a, 6, have = +1,-1, +1/2, —1/2 
respectively and the hamiltonian is block-diagonal with subspaces of fixed total Jz while 
Qi, Q\ change Jz by +1 and Q2, Q\ by —1. 

Consider for instance a supermultiplet containing a = boson \x) of non- vanishing 
mass. Since its p+ and p_ are both non- vanishing this state can be at most annihilated by 
one out of Qi, Q\ and similarly by one out of Q2-, Q\- Acting on \x) in all possible ways we 
see that we generate a supermultiplet containing two bosons with opposite parton number 
parity and two fermions with opposite parton parity. However CPT invariance requires 
the existence of a similar multiplet of antiparticles and we end up with 4 bosons and 4 
fermions as the minimal size of a massive supermultiplell^. 

As an example, consider the anomaly (or Konishi) supermultiplet S = e'^^WaWp. Its 
lowest component is the gaugino bilinear e'^^XaXp which is just proportional to A1A2. Since 
in our setup Ai is expressed in terms of A2 and A we see that, at first order, 5*10) gives a 
state of the type \gga). Such a state is clearly annihilated by Q2 (but not by Q2) and one 
can consistently assume that it is annihilated by Qi (and not by Qi), which mimics exactly 
the situation in D = 4. Acting on \gga) with Qi, Q2 and [Qi,Q2\ (the anticommutator 
being zero), we generate two fermionic states (whose lowest p components look like \ga) 
and a combination of \ggf) and \gah)) and one more scalar (a combination of \ah) and 
\fg) with additional p > 2 components). This gives a total of 2 bosons and 2 fermions to 
which we have to add a similar set of states starting from S\Q) ~ \ ffb). 

As we shall see, our approximation to the full Hamiltonian breaks SUSY to an A/" = 
(1, 1) subgroup of the full J\f = (2, 2) and therefore to the breaking of some degeneracies. 
Nonetheless, an approximate full degeneracy will be seen to remain even in our very simple 
system with Coulomb interactions only. 

^This is only apparently in contrast with the counting of states in 13 = 4 but it is not since, in D = 2 
there is a further doubling of states due to to the distinction between left- and right-movers |19) . 
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3 Cancellation of leading infrared divergences and a 
Coulomb approximation 



In our LC quantization the supercharge Q2 and the momentum P_ are hke in a free 
theory and thus trivial. By contrast, the supercharge Qi and the "Hamiltonian" P+ are 
highly not trivial: non-linear and even non-local. Let us discuss some of their properties 
before making any approximation. From ([9]) we see that every term appearing in Qi is 
trilinear in creation/annihilation operators. Furthermore, since the LC momenta are all 
positive, momentum conservation implies the absence of pure creation or pure annihilation 
terms. Thus Qi connects states with opposite parton-number parity. The LC Hamilto- 
nian, will have either quadratic terms (with one creation and one annihilation 
operator), or quartic terms. Since no pure destruction or pure creation operators appear 
in H^^ the Fock vacuum is annihilated both by the SUSY charges and by the Hamiltonian 
and is an exact zero-energy ground state. The quartic terms induce either 2 — > 2 or 
1 — )■ 3, 3 — )■ 1 transitions and thus conserve (—1)^. 

The full Hamiltonian exhibits, at least superficially, both linear and logarithmic in- 
frared (IR) divergences when some momenta (or combinations thereof) go to zero. The 
logarithmic divergences, that resemble those of the AD theory, presumably need a Block- 
Nordsieck treatment that we plan to implement (analytically and/or numerically) in a 
forthcoming paper. The linear. Coulomb-like divergences are instead neatly cancelled for 
colour-singlet states as we shall now argue. Nevertheless, the finite effects they leave 
behind are expected to dominate the Hamiltonian at large distances. 

In order to illustrate this feature we introduce a rather crude "Coulomb approxima- 
tion", in which we keep only the linearly-divergent terms in the Hamiltonian in their 
minimal form of poles. Such IR singularities only appear in the quadratic terms 
(which are necessarily diagonal) and in the elastic 2 — )■ 2 scattering terms. In the above- 
mentioned approximation the former take the simplified form: 

= -Y.^r Tdk f ^AlA, ; A = /AT, , (12) 
T'" ^ Jo Jo Q 

where A stands for any one of the four parton species. Similarly the elastic quartic terms 
can be simplified as: 



h: 



elastic 
C 



^ 27r 

A,B 



dpidp2 



Jo H 



+ 



P2 



dq 



r 



Ti{AlBlBp,-,A,,^,) 



(13) 



where A, B stand for any one of the four parton species. Note that these elastic terms 
neither change the number of partons nor their species. For this reason we can restrict 
ourselves to subspaces of the entire Hilbert space of states with a given set (i.e. number 
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and species) of partons and fixed total momentum P-{= P below). Furthermore, as 
discussed in Section 1, the 't Hooft limit, N ^ oo with A fixed, selects single trace states 
as the only relevant ones (transitions with multi-trace states being 1/A^-suppressed) . 

In conclusion, we can diagonalize He by splitting the total Hilbert space generated 
by a generic linear superposition of single-trace states into subspaces of definite total 
momentum P and definite parton number p: 

= = Tr |^oI(fci)o^(fc2)...oJ_i(Vi)4(^-E^^)) I0)| ' (14) 

where o\ can be any of our creation operators and all the momenta ki satisfy < ki < P 
for i = I, ...,p and we define kp = P — Yl^=i ^i- Starting with section 5 we will study in 
detail, both analytically and numerically, the eigenvalues and eigenvectors of He restricted 
to ?^p, Hp and Tip. For the numerical approach it will be convenient to discretize the LC 
momenta by compactifying the coordinate x~ . This will be discussed in the next section. 

Before turning to actual calculations let us give a general argument for the cancellation 
of the Coulomb divergences for a general state of the form f|T^ . Let us take an arbitrary 
pair of neighbour partons in f|T^ and consider four distinct contributions to their mutual 
and self-interaction. The self interactions come from the quadratic terms f|T2|) for each one 
of the neighbour partons. However, in order to keep the book-keeping right, we attribute 
half of the quadratic term acting on each parton to its interaction with the left neighbour 
and half to the one with its right neighbour. The book-keeping is easier for the quartic 
Hamiltonian where one just keeps the terms corresponding to the exchange of quanta 
between the two selected partons (there are two such contributions, in general, depending 
on which of the two partons gains energy in the process). 

When these four contributions are added one finds, for each pair of neighbours, the 
following result: 

rA-B,„^ „ „ r^~' ^y{---PA,PB,---)-^{---PA-q^PB + q,---) 



Hc^i- ..PA,PB,...) = — dq- — 



(15) 

where A — B stand for a pair of neighbours and ip{pi, ...,pj\i) is the wavefunction in 
momentum space of the particular chosen state. The total Hamiltonian is given then by 
a sum over A and B. We see that the numerator on the r.h.s. vanishes at g = 0. If we 
take a smooth linear superposition of parton-momentum eigenstates, the numerator will 
vanish quadratically in q thus completely removing the singularity. 

As a physical example of such a smooth superposition let us consider two partons 
displaced by some Aa; = Ax~ in LC space. The momentum wavefunctions will combine 
to give an e**^^^ for the second term in the numerator so that, for each pair: 

A , (l-cos(gAx)) X f°° , sin^{qAx/2) A 



Ht"^^ dq- Y^^^- dq ^1^^ = |Ax| , (16) 

at large separations. We have just discovered that the Coulomb Hamiltonian, acting on a 
single trace state, produces an energy proportional to the sum of the |Ax^ | distances for 
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each pair of neighbouring partons. The constant of proportionahty turns out to be just 
^ i.e. nothing but the 2D string tension for sources in the fundamental representation 
[201 121 j^. This is as it should be, since each parton behaves, vis-a-vis of its neighbours, like 
a parton in the fundamental representation. Its belonging to the adjoint representation 
reveals itself in the fact that it has two neighbours. For a two parton system this just 
produces the well known factor two difference between adjoint and fundamental tension 
in the large- limit. A numerical verification of this analytic argument will be given in 
Sect. 6. 



4 Light-cone compactification 

We now further compactify, for computational convenience, the x~ direction on a circle 
of radius R (with periodic boundary conditions). As a consequence, the p_ momenta are 

h 

P- = n— = ne , n = 0, 1, . . . (17) 
R 



quant ize^^' 



The total conserved momentum is taken to be 

P^ = P = K^ = Ke (18) 
R 

and we shall be interested in the decompactification limit i? — )■ oo (e — )■ 0). In order to 
keep P fixed this also means K ^ oo. Since in LC quantization all momenta are positive, 
P effectively plays the role of an ultraviolet cut off A[/y. In other words, in eq. ffTTI) . 
n < K. If we also take out the zero mode (n = in (ITTI) ). e = h/R plays instead the role 
of an IR cutoff and K = The picture is similar to that of a (2-dimensional) lattice 
gauge theory in Hamiltonian formalism (continuous time and discretized space) with a 
lattice spacing given by a = 27Th/P = L/K and a total of K lattice points to cover a 
large circle of circumference L = 2ttR = Ka. 

According to the LC quantization philosophy, P_|_ plays the role of the Hamiltonian 
and the Lorentz-invariant eigenvalue equation reads: 

M^\s) = 2P+P_\s) = 2PP+\s) = 2KeP+\s) (19) 

The eigenvalues (and eigenvectors) of this operator should have a finite limit a,s K ^ oo. 
We have found it convenient to work in units in which e = 1 (integer parton momenta). 
In that case the operator KPj^ should approach a finite limit as — )■ oo. More physically 
we could have chosen units in which a = 1 and momenta are quantized in units of 2tt /K. 
One can easily check that the limit K oo provides in both cases the same spectrum 
for M2. 



•^This can be generalized to finite N with Cp = {N"^ — 1)/2N replacing N/2 (private communication 
by G. C. Rossi). 

"^In this section we shall keep Planck's constant explicit in order to illustrate the emergence of a 
quantum-mechanical string. 
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Figure 1: Cutoff dependence of the eigenvalues of He in the three lowest multiplicity 
sectors. The size of the He matrix grows with the resolution parameter K and can be 
read from the range covered on the horizontal axis, K = 30, 50, 70 for p = 2,3, and 
K = 15, 25, 35 for p = A. 

We expect to find finite-size effects at finite K. Because of periodicity the maximal x~ 
distance is Ka/2 and therefore only states whose wavefunction is concentrated in regions 
much smaller than Ka/2 are expected to have already reached an asymptotic limit at 
some given K. We will see later that this is exactly what happens numerically, but let 
us anticipate the presentation of some of the results we will find for p = 2. Once we have 
discretized the momenta in units of e = 1, the operator becomes an [K — 1) x [K — 1) 
matrix whose low eigenvalues converge fast to some finite values while increasing K. This 
can be seen in Fig. [1] (to be discussed in detail in Section [6]) which illustrates what 
happens for K = 30, 50 and 70. This is equivalent to discretizing states with a fixed 
total momentum P = 1 by subsequently taking e = 1/30, 1/50 and 1/70, We see a neat 
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convergence of the lowest eigenvalues towards a smooth (and almost linear) spectrum. 

Referring still to the two-parton sector, we shall also find that the eigenvectors exhibit 
a sharp localization in position (i.e. x~) space, or, more precisely, in the distance Ax~ 
between the two partons. The average distance appears to be quantized: 

{Ax-)n = Cn^ = 2TTCn^ , n = 1, 2, . . . - 1 , (20) 

where the c„ are a sequence of numbers going from a number of 0(1) to a number of 
0{K/2) in steps of 0(1). The ordering coincides with that of the energy eigenvalues. In 
other words, the average distance is of the order of a = 27ih/P_ for the lightest states and 
grows up to the maximal allowed physical distance for the heaviest eigenstates. As we 
increase K, the low-lying states stabilize while the heavy ones change and stabilize only 
at higher values of K (i.e. when Aa;~ is well within the compact circle). The spread in 
Ax~ is 0(a) for all the states (this is the above-mentioned position-space localization). 

The low- lying eigenvalues of P+ behave like Pl^ for large P_ and, indeed, one finds 
(c.f. Figj9]) an approximate linear relation between energy and average distance: 

Pj") ^ A(Ax-)„ (21) 

where A is 't Hooft's coupling normalized in such a way that it corresponds to a classical 
tension. This means that an appealing string picture emerges whereby energies are pro- 
portional to the string tension A. This result, however, is Lorentz- frame dependent. In 
order to find a Lorentz invariant result we compute the mass eigenvalues: 

= 2P_P|"^ ~ 2AP_(Ax-)„ = AnhXcn = AM„(Aa;""')„ (22) 

where we have used first fl20|) and, for the last step, an ra-dependent Lorentz transformation 
with boost a/2P_/M„ in order to go to the nth state rest frame. Thus we finally obtain: 

(Ax^™)„ = Vi^/, , M„ ^ A(Ax^"^)„ = Vi^M, (23) 

in terms of the string length and mass scales: Is = \/h/X, Mg = \/hX. This is in 
perfect agreement with expectations |20[ 121] once we realize that, for p = 2, the sum over 
neighbour-parton distances is 2Ax. 



5 Bases and matrix elements at large Nc 

Below we quote a few explicit expressions for matrix elements of the discretized (and 
rescaled) Hamiltonian he defined by: 

He = -^hc. (24) 
The mode expansion of he in terms of the discretized creation and annihilation operators 

= f-pr^me (25) 

V R 
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is the same as in fll3|12p with all momenta p ^ m and integrals J dp ^ Ylm- 

There are four kinds (species) of partons in our model: two bosons and two fermions 
with the corresponding annihilation operators denoted by a,b,f,g. The detailed con- 
struction of the planar bases in each multiplicity sector differs slightly depending on the 
particular choice of parton species involved and the same applies to the matrix elements 
of he ([TTO]). 

5.1 Two different partons 

At given resolution K the states belonging to the discretised version of Tij^ are labeled 
by one integer: 

\n) = Tr[AKB^K-nm, n=l,...,K-l, (26) 
where A,B = a, b, f, g. The large N rules, developed e.g. in [22], give 

{n\hc\n') =6n,n' {C{n) + C{K - n)) - (1 - 5„,„0 ^ 



{I - 5K~n,K~n')—, , (27) 



n—l ^ 



5=1 ^ 

5.2 Two identical partons 

In this case, due to the cyclic symmetry n K — n, only half of the states from above, is 
linearly independent. The matrix elements can be neatly expressed by ( ant i) symmetrized 
ones from the previous case 

K -1 

{n\hc\n')id = {n\hc\n')dif ± {K - n\hc\n')dif , n=l,... — - — , K odd, (28) 

where the suffix id stands for identical partons and the suffix dif for different partons while 
the + sign is for identical bosons and the — for fermions. The union of the two (bosonic and 
fermionic) spectra for identical partons will reconstruct the complete spectrum obtained 
from eq.([2j 



5.3 Three different partons 

We consider now states composed by three different partons (taken as abf). These are 
labeled by two integers: 



n, m 



)=Tr[alblfj,_^_J\0), l<n<K-2, l<m<K~n-l (29) 
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and the matrix elements of Eq. fll2lll3p read: 



(m, n\hc\Tn\ n') =5rn,rn' ^n,n' {C {m) + C{n) + C{K — m — n)) 

- - 5m,m')/2("^ - m'f 

- 5m+n,m'+n'(l " '^m,m')/2("^ - m'Y . (30) 

In the case of states with just two parton species (e.g. aaf) he has the same matrix 
elements. The Hilbert space is however twice as small since cyclic and anticyclic permu- 
tations correspond to the same state. On the other hand the bigger Hilbert space for 
three different species splits into two separate sectors, since states corresponding to above 
cyclic and anticyclic permutation do not interact in the planar limit. 

5.4 Three identical partons 

With a trick similar to the two parton case one can exploit the cyclic symmetry and 
express matrix elements in terms of those for three different partons 

(n, m\hc\n\ m')id ={n, m\hc\n , m')dif + {m, K - n - m\hc\n , m)dif+ 

+ {K - n - Tn,n\hc\n,m)dif ■ (31) 

This time there is no difference between bosonic and fermionic sectors since Z3 shifts 
involve an even number of fermionic transpositions. 

5.5 Four different partons, the abfg sector 

The planar basis is now parametrized by three integers 

|n, m, o) = Tr[aJ,6|„/]c/^_„_„_J |0) , 

l<n< K -3, l<m<K-n-2, 1 < o < K - m - n - 1 , 
and matrix elements of Eq. fll2lll3p read 

(m, n, o\hc\m' , n', o') =Sm,,m'Sn,n'So,o'{C{^) + C{n) + C(o) + C{K — m — n — o)) 

— 5m,m'5n,n'(l ~ ^o,o')/2(o " o'Y 

- ^n,n"5o,o'(l - 5m,m')/2("^ " ni'Y 

~ ^o,o'^m+n+o,m'+n'+o'i^- " ^n,n') f^i^ — Tl)'^ 

~ Sm+n+o,m'+n'+o'Sm,m'{^ ~ '^o,o')/2(o — O )^ . (32) 

Generalization to higher parton multiplicities is pretty straightforward. 
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6 Results 



We turn now to discuss quantitative solutions obtained, mainly numerically, with the 
aid of Mathematica. The crude Coulomb approximation introduced above is not only 
the lowest approximation for the full set of the complete coupled LC eigenequations but, 
at the same time, it defines a natural generalization of 't Hooft equation to many-body 
sectors. Rather than concentrate on separate parton multiplicities, p, we shall focus on 
a few physical issues and compare results for different multiplicities. Until now only the 
first three nontrivial sectors, i.e. p = 2,3,4 have been looked upon. Obviously, further 
increase of p is technically more challenging, however it is feasible, if such a need arises, 
by employing more dedicated methods and algorithms. 

6.1 K dependence vs. K = oo limit 

Figured] summarizes the cutoff dependence of spectra of He for different multiplicities. As 
expected, lower states converge faster, however one should remember that at the high end 
of the spectrum new states appear for each K. Since K controls the length L = 2ttK/P 
of our torus, the highest state, for example, will never converge because it is a new state 
with higher and higher energy as we increase K. 

The eigenenergies of PHq (that is are displayed as a function of a single index n 

which labels consecutively ordered eigenvalues. Only for the two parton case n is directly 
related with a single physical observable (see below). For more partons, n is effectively 
composed of more quantum numbers related to other, yet unknown, quantities conserved 
in a particular many body sector. As a consequence, the p = 2 dependence is nicely linear, 
at large K, while the limiting curves for higher p are more complicated. They reflect the 
above degeneracy with more and more states below a fixed energy as we increase p. In 
fact one can read the spectral density of states dN/ dE directly from the figure. It grows 
with E at intermediate energies, the growth becoming more and more rapid as we increase 
p. However at the highest end of the spectrum the density saturates revealing some sort 
of blocking related to periodicity. 

Let us now look for the cutoff dependence of the lowest levels in more detail. Figure [2] 
shows the first few eigenenergies in the boson-boson and fermion-fermion, p = 2, sectors as 
a function of 1/K. Evidently dependence on K is very weak and the extrapolation to the 
continuum momentum limit is straightforward. In fact, in this limit the LC eigenequations 
are nothing but 't Hooft equations with adjoint charges. Hence the limiting values can 
be obtained by solving directly these equations. They are also displayed in Figure [21 
providing a rather satisfactory cross check of the whole approach. 

Finally we comment on the difference between bound states made of fermionic and/or 
bosonic identical partons. Since He is invariant under the reflection m -ir^ K — m, the 
non-degenerate eigenstates have definite parity. This is the case in the two parton sector. 
Moreover, the Z2 symmetry of planar states together with the "flavour" independence al- 
lows to identify the even and odd solutions as identical boson-boson and identical fermion- 
fermion bound states. In the K = 00 case, we have generated corresponding eigenstates 
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Figure 2: Convergence of the lowest energy levels with the cutoff (100 < K < 3000) 
and comparison between the even (identical boson-boson) and odd (identical fermion- 
fermion) solutions. Stars represent direct solutions of the 't Hooft equations for continuous 
momenta [K = oo). 

by employing bases with the required symmetry. 

As usual, the lowest state is symmetric. It's wave function is constant in the parton 
momentum, k, with exactly zero eigenvalue for any cutoff K^. This state is not displayed 
in the right panel of Figure [2l 

Let us quote, for completeness, the continuum momentum limit of the eigen-equations 
of He which was used to obtain the 1/K = points in Fig. [21 They can be readily 
derived from the discretized matrix elements fl27ti32p or, equivalently, by applying our He 
f lT2|T3|) to the n-parton state and using rules of the planar calculus [22]. For two partons 

^Many other massless states have been found in the hterature 18 for any K . These are due to a 
finite- artifact by which states with a large number of partons, k ^ K, are annihilated by the SUSY 
charges (and therefore by the Hamiltonian) , since all transitions are blocked by the finite momentum 
resolution. Obviously only states with fc <C if can be reliably computed in our approach. 
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one obtains: 



Ec^{p) 



2^ 



(33) 



which indeed is equivalent to 't Hooft equation: 



EMx) = ^PV f dy 



(X - y)2 



(34) 



It is perhaps worth observing that in the former variables, i.e. in terms of the momentum 
transfer g, the linear and logarithmic divergencies clearly cancel leaving behind the finite 
right hand side of ( l33l) . This provides the justification for the principal value prescription 
commonly used in (!34|) . 

The multiparton generalizations follow from ( !T2|T3|) as well 



E^{pi,...Pn) = 

^ n p^_, 99(pi,...j9j,pj+i,...p^) -(p{pi,...pi -q,Pi+l +q,...pn) 

i = l ^J-P^ + l+^ y 



(35) 



and, again, can be written in many equivalent forms [23] . 



6.2 Spatial structure of multiparton states 

It turns out that the eigenstates of Hq have a very simple and natural structure which 
shows up most beautifully in position space. For p = 2 it is summarized in Fig. [3] where a 
sample of two parton eigenstates spanning the whole interval of eigenenergies is displayed. 
What is actually shown is the modulus squared of the discretized version of the Fourier 
transform 

V'.(Ai2)= r e-'^''P^MPi,P-Pi)dpi, (36) 
Jo 

as a function of the relative, light cone, distance A12 = Xi — X2, with ipr being the 
r-th eigenstate of ( 1271) . Upon the discretization all momenta and coordinates become 
dimensionless integers: P ^ K, pi ^ m, A12 — )■ du- 

The interpretation of the above Figure is clear. The eigenstates are very well localized 
in relative distance, and there is a very strong (linear in fact) correlation between du and 
r. In the lowest state two partons sit on top of each other and their energy is exactly 
zero. Excited states correspond to partons separated by a finite distance which gradually 
increases with the energy. Finally, in the highest state partons are maximally separated, 
i.e. are located at the antipodes of the circle (the Figure refers to A' = 201). The relation 
between di2 and eigenenergy turns out to be linear, as expected in the one dimensional 
world (c.f. Sect. 6. 4). This also explains the linearity of the p = 2 plot in Fig. HJ 
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Figure 3: Space profiles of various eigenstates ( l36l) witli two partons, K=201 (r = 
1,26,50,100,150,200). 

Interestingly, this picture generalizes naturally to many parton sectors. Figure H] dis- 
plays a sweep through three parton eigenstates obtained by diagonalization of ( 130|) . Co- 
ordinate space wave functions depend now on two independent relative distances which 
we choose as A13 = xi — X3 and A23 = X2 — x^,. 

7/>,(Ai3, A23) = / e'^»^'^e^^^^^^V>,(pi,p2, P-Pi- P2)dpidp2 . (37) 

Again, after the discretization the relative distances become integer: 1 < di^, ^23 < K — 2, 
with K — 2 being the period of the discrete Fourier transforms. 

Similarly to the p = 2 case the wave functions are composed of a series of very narrow 
(in lattice units) structures which give a sharp localization in relative distances. Again 
the energy of the lowest state is exactly zero with all three partons located at the same 
point (upper left panel). Going to higher and higher states partons are moving apart 
migrating into the whole circle and finally, in the highest state, three partons are sitting 
at the maximal and equal distances forming the familiar "mercedes" star (lower right). 

Obviously the structure of three parton states is much richer than that in the p = 2 
sector. Nevertheless a number of regularities can be found which provide a compelling 
overall picture. They are better seen in the contour plots which we will now discuss. 

Apparently the whole spectrum consists of a series of states with similar properties. 
Beginning of one such series is shown in FigjSl Again we see that the higher the energy, the 



16 



larger are the inter-parton distances. The new feature is that now there are many peaks 
(i.e. parton configurations) in a single state. Part of it comes from the Z3 symmetry 
of the displayed densities. However the rest provides a beautiful confirmation of the 
linearity of the Coulomb potential and/or the underlying string picture. Namely, all 
parton configurations, composing a particular state in one series, appear to have the 
same value of the "combined string length" I = \di2\ + \d23\ + Ic^sil- And vice versa: 
configurations with different / belong to states with different energy. As an illustration 
consider the fifth state displayed in Fig. [51 It has / = 12, this can be achieved, e.g. 
with partons 1 and 2 separated by 6 units and with parton no. 3 somewhere between the 
two. This is represented by 5 configurations (peaks) extending from {dis,d23) = (1,-5) 
to (5, —1) as no. 3 moves from no. 1 to no. 2. Similarly one can decode all structures 
appearing in that Figure. One more example: a horizontal ridge extending between 
(—5, —6) and (—1, —6), in that panel, describes partons 2 and 3 separated by 6 units and 
parton 1 located in five positions between the two. 




Figure 4: As in Fig|3l but for three partons and in terms of two relative distances ^13 and 
d23, r = 1,80,120,406, K = 30. 
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Figure 5: A clean series (A) of three parton eigenstates, K = 100, 4 < Z < 14. 

The Z3 symmetry of planar states is not evident in these figures, because of the 
asymmetric variables used. However it is there and corresponds to {x,y) = (^13,(^23) — )• 
{y — x,—x) — ?► {—y,x — y). The Z3 symmetric representation is shown in FigE] where 
another series is displayed on the massless "Dalitz plot" ensuring the constraint di2 + 
^23 + dsi = 0. 

Finally, let us comment on the existence of different series in the many body spectrum 
which we find quite intriguing. As already said, we see "experimentally" that our three 
parton states group naturally into series. States in one series differ only by the increase 
of the relative distance between partons, as in Figj5l However states from various series 
exhibit other differences. An example of another series was just given in Figj6l Here again, 
increasing the inter parton distances increases the eigenenergy, however the pattern of the 
configuration remains the same. On the other hand, apart from the change of a display, 
there is a clear difference between patterns shown in Figures [5] and [61 In the first series (A) 
the partons never coincide, while in the other one (D) in every configuration one relative 
distance vanishes. However, we have also found other series where the differences are not 
so clear. In general, states from different series can have the same combined string length 
/ and yet they have different energies. This suggests the existence of other conserved 
quantities, hence quantum numbers, which also control multiparton spectra. 

Solutions in the four parton sector show qualitatively the same phenomena. In Fig[7]we 
display contour plots, in three relative distances d^, 624, and 6/34, of various eigenstates of 
(!32|) . Again, in the lowest state all relative distances vanish; then, in higher states, partons 
gradually separate and in the highest state they group in two closely bound pairs sitting 
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Figure 6: Different series (D) and on the "Dalitz plot". 

at the antipodal locations on the circle. This is different than the naive extrapolation from 
the three parton case and illustrates how rich is the system with possibly more surprises 
at higher multiplicities. Similarly to p = 2, 3, states are sharply localized in the relative 
distances and separations between various partons can be determined. 

Obviously, multidimensional representations like FigJT] become unpractical for higher 
multiplicities. An alternative way to proceed is to study inclusive densities and correla- 
tions. In our case, e.g. in a given p parton sector, an inclusive single parton density can 
be defined as 

/p-i 
dP-'A,Y,^{A - A„,)|V^,(4)p (38) 
i=i 

and gives the number of partons at a distance A from i.e. the last one. It can be easily 
calculated from our exclusive wave functions or, yet simpler, directly from the Fourier 
components. The latter representation reads, e.g. in the four parton case, 

Dr{A)= f dp2dp3\M^,P2,P3)\^ + cycl. (39) 

Jp2,P3,P-P2-P3>0 

with '?/'r(A,p25P3) standing for the partial Fourier transform - only in the first variable. 

In Fig. [S] above density is shown in the four parton sector. It confirms what we have 
already learned from the exclusive data. However the structure of four parton states 
can now be seen in more detail. Rather than sweeping through the whole spectrum, we 
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Figure 7: Structure of eigenstates with four partons. Contour plots in three relative 
distances (d^, 6/24, ^34) for states no. 1,9,35,60,100,165 spanning the whole range of states 

for K = 12, Vmax = 165. 

have concentrated on lower states. FigJH] clearly shows the growth of the distance be- 
tween the two outermost partons and how the intermediate positions between the two are 
populated. Of course the complete information could be recovered only upon examining 
simultaneously higher inclusive densities. 

6.3 Some analytic considerations 

Let us now discuss some analytic aspects of the solutions. To this date such solutions 
are not available in spite of many attempts [24J. Generalization for many bodies only 
increases the challenge. We believe however that the numerical studies presented here 
can also contribute to the analytic understanding of these systems. Consider for example 
the massless, Ec = 0, bound state. It has been found numerically, and it is obvious from 
( 135|1 that such a solution exists for all multiplicities and has a constant wave function in the 
momentum representation. It is then a simple matter to construct its (LC) configuration 
space counterpart. For two partons we have 



ip{xi,X2) = c / dpidp2e 
Jo 

^ g-jP(x-i+X2)/2^ 





sin(PAi2) 



(40) 
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with the (not normahzed) probabihty distribution depending only on A12 



Kf^5i!^. (41) 

For three partons one obtains analogously 

^(Xi, X2, X3) = C / dpidp2dp3 Sipi +P2+P3- P)e-'Pi"i-'f2"^-'''3"3 = 

Jo 

(Aie-*^^^ sin(A2P) - Aac"*^^^ sin(AiP)) , (42) 



2AiA2(Ai - A2 

where we used Ai = Ai3/2,A2 = A23/2 as relative distances. Upon normalization 
c = P/ v^, we obtain for the density 

+4(Ai siii2(A2F) - A2 sm2(AiF))2] , (43) 

which is completely symmetric in exchanging Ai with A2 and regular for Ai, A2 — ?■ 0. As 
for the two parton case we see that in the limit P — )■ cxd we find a (5-function centered at 
zero, both for Ai and A2, confirming that actually the zero energy state corresponds to 
three partons sitting all together at the same position. 

Generalizing this computation to higher sectors is trivial, hence the massless eigen- 
states can be constructed analytically in all multiplicity sectors. 



6.4 String picture 

The findings reported so far support the widely accepted string picture of (1+1)- dimen- 
sional planar gauge theories. The relation is quite natural, even hardly surprising, in the 
case of two partons. It is also generally expected with more partons, however our results 
provide a direct illustration of how this is happening. 

We have seen that the two parton eigenstates are very well localized in (LC) configu- 
ration space. Consequently, the length of the effective string between two partons can be 
readily extracted from our eigenstates (c.f. FigJSj). Figure [9] shows the dependence of the 
eigenenergies on that length. It approaches a nice linear form at large cutoff, K, with a 
well defined, finite, string tension. 

There is an interesting correlation between the parity under the reflection k i-^ P — k, 
and the inter parton distances. Namely, in even (odd) states partons are separated by 
integer (half-integer) distances. This explains why the eigenenergies of the symmetric and 
antisymmetric states are half a way between each other in Fig. [2] 

Before moving to higher multiplicities, let us compare our numerical results with the 
theoretical prediction for the effective string tension in the two-parton case, a = A/2 
[20l[2T]. FigJTOl shows the ratio of two parton eigenenergies to the combined string length 



22 



100 



A 



40 - 



20 



50 



100 

P(kl2l+U-2ll) 



150 



200 



Figure 9: Two-parton energies, as a function of (doubled) parton separation, for increasing 
K. Different colors (bottom to top) correspond to = 50, 100,200,400. 
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Figure 10: Left: as in Figj9]but for the string tension together with the extrapolations to 
K = oo (black). Right: polynomial extrapolations (see the text) to a = 0. 



X = \xi2\ + |2;2i|, in units of A, as a function of the dimensionless lattice distance Px = 
2TTx/a = 2tiI. Different colors show results for different cutoffs, K = 50, 100, 200, 400, 800. 

As expected, there is a significant cutoff dependence for larger parton separations. The 
K dependence is however rather weak and can be easily taken care of. Black points show 
results of the polynomial (in l/K) extrapolation to i^' = oo, at few values of I. In lattice 
terminology they correspond to the string tension determined from finite lattice distances 
hence are still biased by finite a effects. To get rid of the latter one has to perform the 
continuum (i.e. a — )■ 0) limit. This is summarized in FigJTO] where the a/x dependence of 
above K extrapolations is displayed. Again the a/x dependence is mild and polynomial 
fits provide quite stable extrapolations to a = 0. Those are in very satisfactory agreement 
with the theory cf. Table I. 

With three partons, situation is yet more interesting. Assigning automatically the 
combined string length / to each state on the basis of integer coordinates of sharp peaks 
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M 1 2 3 

WAf(O) 0.4944 0.4980 0.4997 



Table 1: First three polynomial extrapolations of the a dependence of the string tension 
from FigJTOl 
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Figure 11: Eigenenergies of three parton states vs. the combined string length /. 



in the density, results in Fig. [TT](left, K = 100). This only confirms the existence of some 
ambiguities. We have already found, however, that there are series of states composed 
of similar patterns of configurations. Could they account for what is seen in Fig. 1111 .'' 
In fact yes, in the right panel of Fig. [11] we show similar plot but constructed for five 
series identified by inspection of lowest 20 eigenstates. Two of them were described in 
the previous subsection. Indeed, for each series we observe a clean linear growth of the 
energy with /. The slope of that dependence is consistent (albeit somewhat smaller) with 
the one seen for two partons. 
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Figure 12: As in FigJTT]but for four partons. 
Similar behaviour is seen in the four parton sector. Fig. [121 obviously the situation is 
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more complex with more families and larger spread of energies at fixed "1". No attempt 
was made to identify separate series yet and to extract the string tension. The latter 
requires the former and extrapolations in K and a, as was done for p = 2. 

Clearly more detailed studies are needed to unravel a complete structure of higher 
parton states. In particular comparison with the direct solutions of 't Hooft equations in 
higher-p sectors would be very useful. 

6.5 Pre-SUSY 

As already mentioned our drastic approximation breaks half of the original supersymmetry 
leaving behind just the ^{1, 1) subgroup generated by Q2-,Q\- In spite of this it is amusing 
to compare the spectra of sectors that should be connected by the action of all the SUSY 
charges. The supersymmetry generated by Q25 Q2 connects states with the same number 
of partons and is self-evident in our graphs. On top, there are further degeneracies between 
states with the same p due to our Coulomb approximation. It is clear, instead, that there 
is no exact degeneracy for states of different p even if some of them should be connected 
through the action of the Qi, Qi generators. 

aaa aa xxy/xyz ff fff 
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Figure 13: The lowest levels with 2 and 3 partons. Here a represents either a or 6 and / 
either / or g. For p = 2 the aa and // spectra, for identical partons, are complementary 
and their sum gives all the unequal-parton spectra. For p = 3 the xxx (i.e. aaa or ///) 
spectra are all identical but contain far less states than the xxy or xyz spectra. The xyz 
spectra are equal to the xxy spectra modulo a doubling of the states since xyz and xzy 
are two distinct, but degenerate, states. 

Let us consider, as an example, the Konishi (anomaly) supermultiplet keeping Fig. [13] 
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in mind. As discussed in section 2.1 the full chiral supcmultiplet in this case contains p = 2 
states of non-identical partons and p = 3 states in which there are at least two different 
species. The corresponding energy levels match the union of the aa and // spectra in the 
figure and the central (xxy/xyz) levels. We see that there is (decreasingly) good matching 
for the first three excited levels while for the fourth the matching is already quite poor 
(for higher levels the good matching looks a bit accidental given the high density of p = 3 
levels) . 

Consider now a different supermultiplet, one that contains the state \aa). Its partners 
should be found in the af sector and in the aaf and aff sectors and all works as in the 
previous case. 

Consider finally the supermultiplet containing the state |//). Apparently we find 

a problem since the lowest excited // state has no nearly-degenerate partner in the 
/// sector. However this is as it should be: the \ ff) wavefunction must be odd under 
the interchange of the two momenta. When one applies Qi to it one finds that this 
antisymmetry clashes with the symmetry needed in \ fff). Instead, the \ ff) state should 
form a multiplet with \af), \afg) and \abf). The matching is now excellent. 

In conclusion, although a priori our drastic truncation of the Hamiltonian could have 
left no sign at all of the supersymmetry generated by Qi and Qi, we have found that 
some trace of the full J\f = (2, 2) supersymmetry appears to have survived in the spectra 
of He- This makes us confident that our truncated Hamiltonian may represent a fairly 
good approximation to the exact one. 

7 Conclusions and outlook 

We have considered the dimensional reduction of D = 4, A/" = 1 SYM theory to two 
dimensions in the large-A^ (planar) limit using the LC gauge and LC quantization. This 
allows us to explicitly eliminate all non-physical degrees of freedom at the price of having 
a non-local LC Hamiltonian. Nonetheless, such an Hamiltonian exhibits many desirable 
features: it is invariant under N' = (2, 2) supersymmetry transformations provided these 
are suitably defined in order not to destroy the LC-gauge choice; it is manifestly normal- 
ordered and positive-semidefinite so that it possesses an exact zero-energy ground state 
and a spectrum of non-negative energy excitations. 

In order to solve for the eigenvalues and eigenstate of the Hamiltonian we have com- 
pactified LC-space to a circle of radius R allowing to work within a finite Hilbert space 
as long as we keep the conserved LC momentum of our states finite, the idea being, of 
course, to eventually send R to infinity and check that physical quantities approach a 
finite smooth limit. 

Infrared (IR) divergences appear to make this task somewhat technically complicated 
(although in principle possible) and therefore, in this first paper, we have truncated 
the Hamiltonian to what looks superficially as its most IR divergent part. Indeed, for 
the colour-singlet, single-trace states that survive in the large- limit, these linear IR 
divergences are neatly cancelled and get replaced by an effective Coulomb interaction 
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which, in D = 2 gives a confining hnear potentiaL We have then studied many properties 
of the eigenvalues and eigenfunctions of this Coulomb Hamiltonian, Hq, and confirmed 
that, at least in this approximation, the large R limit is smooth (although larger and 
larger R are needed before heavier and heavier states stabilize). 

The resulting model breaks supersymmetry down to an A/" = (1, 1) subgroup and looks 
like a supersymmetric generalization of 't Hooft's original model [12] where states with 
an arbitrary number of partons are present even at leading order in 

When seen in position-space a nice string picture emerges in which the mass of each 
state is proportional to the sum of the (center-of-mass distances) of each pair of neighbour- 
ing partons. These distances are quite sharply quantized leading to a discrete spectrum 
with approximately linear "Regge" trajectories. The numerical value of the proportion- 
ality constant (the string tension) agrees very well with theoretical expectations. For 
states with more than two partons the situation is obviously richer. Yet we were able to 
identify clean series of three-parton states whose eigenenergies are indeed proportional to 
the combined length of strings stretched between neighbouring partons. Moreover, the 
detailed patterns of parton configurations contributing to these states confirm unambigu- 
ously the linear form of the two-body interactions. String tension seen in these sectors is 
compatible with the one extracted from the two parton sector. 

In our Coulomb approximation the first excitation over the Fock vacuum is also mass- 
less with all the partons sitting at the same point but this is most likely an artifact of our 
Coulomb approximation that allows all partons to sit at the same point without paying 
any kinetic-energy price. Since even those massless states are nicely paired in supermul- 
tiplets, we expect them to be lifted to some finite energy. Indeed, if we add some finite 
terms present in the full Hamiltonian (|5]), we see that these state acquire a non-zero en- 
ergy. At least in our approximation we find no evidence (apart from the just mentioned 
states) for the absence of a mass gap reported in some previous studies [IS]. We also see 
clearly how the existence of many other massless states, for each value of R, is nothing but 
a consequence of the breakdown of the method when the number of partons approaches 
its maximal value compatible with momentum conservation. 

Finally, we found that, even in our drastic approximation the full M = (2,2) super- 
symmetry of the original model shows up as an approximate supermultiplet structure at 
least for the lightest states. Our belief is that the "worst" divergent part in gives us the 
bulk structure of the energy states, namely the nice discrete linear spectrum. We expect 
the logarithmic IR singularities to lead to a dressing of our states a la Block-Nordsieck 
and to modify the energies by lifting, in particular, the zero-energy states. We plan to 
report on progress in this direction in a forthcoming paper. 
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Note Added 

After having circulated a preliminary version of this paper we became aware of previous 
work [25j, [22], [27], [2E], claiming that two-dimensional theories of the kind we considered 
(i.e. with massless fermions in the adjoint representation) should exhibit a Schwinger-like 
phenomenon even in the planar large- limit. This would result in the screening of the 
linear potential and in the spectrum becoming continuous above a certain energy scale. If 
this were true, it would mean that the terms we neglected should have a dramatic effect, 
at least for the high-energy part of the spectrum. This point, as well as the dependence 
of the spectrum itself on the scalar VEV, clearly deserve further investigations. We are 
grateful to A. Armoni for bringing the abovementioned papers to our attention. 
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